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GENERATION OF A DESIRED THREE-DIMENSIONAL ELECTROMAGNETIC 
FIELD 

FIELD OF THE INVENTION 

The present Invention relates to a method and a system for synthesizing a prescribed 
5 three-dimensional electromagnetic field based on generalized phase contrast 
imaging. 

BACKGROUND OF THE INVENTION 

It is well known to form an Image by phase contrast imaging methods in which phase 
modulation of light is converted into intensity modulation. As opposed to intensity 
10 modulation, phase modulation does not involve loss of energy. 

A generalized phase contrast imaging method and system for synthesizing a 
prescribed intensity pattern is disclosed In WO 96/34207, which is hereby 
Incorporated by reference. The generalized method is not based on the so-called 
Zemike approximation that the phase shift <p is less than 1 radian. An improved 
15 method is provided without this approximation and based on imaging with a simple 
one-to-one mapping of resolution elements or pixels of a spatial phase modulator 
and resolution elements of the generated intensity pattern. 

The disclosed phase contrast imaging method of synthesizing a prescribed intensity 
pattern /(tfy), comprises the steps of: 

20 • dividing the intensity pattern I(x\y>) =\o(x\ytf into pixels in accordance with 
the disposition of resolution elements (x,y) of a spatial phase mask having 

a plurality of individual resolution elements (x,y), each resolution element (*,y) 
modulating the phase of electromagnetic radiation incident upon it with a 
predetermined phasor value e lmi t 

25 • radiating electromagnetic radiation towards the spatial phase mask, 

• Fourier or Fresnei transforming the modulated electromagnetic radiation, 

• phase shifting with a spatial phase filter (SPF) In a region of spatial frequencies 
comprising DC in the Fourier or Fresnei plane, the modulated electromagnetic 
radiation by a predetermined phase shift value 9 in relation to the remaining part 

30 of the electromagnetic radiation, and 

• forming the intensity pattern by Fourier or Fresnei transforming, respectively, the 
phase shifted Fourier or Fresnei transformed modulated electromagnetic 
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radiation, whereby each resolution element (x,y) of the phase mask is imaged on 
a corresponding resolution element (^.y 1 ) of the image, 

• calculating the phasor values e^ of the phase mask and the phase shift value 
0in accordance with 

for selected phase shift values 0, a being the average of the phasors rf** of 
the resolution elements of the phase mask, 

• selecting, for each resolution element, one of two phasor values which represent 
a particular grey level, and 

1 0 • supplying the selected phasor values rf** to the resolution elements (x t y) of the 
spatial phase mask. 

In one embodiment disclosed in WO 96/34207, the spatial phase mask is positioned 
at the front focal plane of a lens white the spatial phase filter is positioned in the back 
focal plane of the lens, whereby a first electromagnetic field at the spatial phase 

15 mask is Fourier transformed by the lens into a second electromagnetic field at the 
spatial phase filter Thus, specific spatial frequencies of the first electromagnetic field 
will be transmitted through the spatial phase filter at specific positions of the phase 
filter. For example, the energy of the electromagnetic radiation at zero frequency 
(DC) Is modified and transformed onto the intersecting point of the Fourier plane and 

20 the optical axis of the lens also denoted the zero-order diffraction region by the phase 
filter. 

SUMMARY OF THE INVENTION 

The above-mentioned method operates on a plane Incoming electromagnetic field 
with the aim to generate a two-dimensional intensity pattern. It is an object of the 
25 present invention to provide a method and a system synthesizing a prescribed three- 
dimensional electromagnetic field for further flexibility, for example to be able to focus 
light energy in a three-dimensional volume. Such a method and apparatus may be 
utilized in advanced optical micro- and nano-manipulation, such as by provision of a 
multiple-beam optical tweezer. 

30 According to a first aspect of the invention the above-mentioned and other objects 
are fulfilled by a phase contrast system for synthesizing an output electromagnetic 
field u^', y'\z"), comprising 
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a first phase modifying element for phase modulation of an input electromagnetic 
field by phasor values 

first Fourier or Fresnel optics for Fourier or Fresnel transforming the phase 
modulated electromagnetic field positioned in the propagation path of the phase 
5 modulated field, 

a spatial filter for filtering the Fourier or Fresnel transformed electromagnetic 
radiation by 

in a region of spatial frequencies comprising DC in the Fourier or Fresnel 
plane 

1 0 phase shifting with a predetermined phase shift value 9 the modulated 

electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constants, and 

15 in a region of remaining spatial frequencies In the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant A, 

second Fourier or Fresnel optics for forming an electromagnetic field o(x\ /) by 
Fourier or Fresnel transforming the phase shifted Fourier or Fresnel transformed 
20 electromagnetic field, and 

a second phase modifying element for phase modulating the electromagnetic field 
o(x\ y 1 ) into the output electromagnetic field o(x\ M propagating as the 
desired output electromagnetic field u(x\y\z n ). 

According to a second aspect of the present invention, the above and other objects 
25 are fulfilled by a method of synthesizing an output electromagnetic field u(x", y", z"), 
comprising the steps of 

phase modulating an input electromagnetic field by phasor values e im \ 

Fourier or Fresnel transforming the phase modulated electromagnetic field, 

filtering the Fourier or Fresnel transformed electromagnetic radiation by 

30 in a region of spatial frequencies comprising DC In the Fourier or Fresnel 

plane 
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phase shifting with a predetermined phase shift value 9 the modulated 
electromagnetic radiation In relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
5 with a constant B, and 

in a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant A, 

forming an electromagnetic field o(x\ y') by Fourier or Fresnel transforming the phase 
1 0 shifted Fourier or Fresnel transformed electromagnetic field, and 

phase modulating the electromagnetic field o(x\ y 1 ) into the output electromagnetic 
field o(x\ f)Q w V.f) propagating as the desired output electromagnetic field 

The method may further comprise the steps of 

1 5 dividing the electromagnetic field ofx'.y 1 ) into pixels in accordance with the disposition 
of resolution elements (x, y) of a first phase modifying element having 

a plurality of individual resolution elements (x, y), each resolution element 
(x, y) modulating the phase of electromagnetic radiation incident upon it with 
a predetermined phasor value e^, 

20 calculating the phasor values of the phase modifying element and the phase 
shift value 0 substantially in accordance with 

o{x\y) Si*[exp(^(x , B y , ))+ JTpKw-' exp(/0) -l)] 

wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
25 order diffraction region, 

B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region, 

a a [ff|exp(/^) is the average of the phasors of the resolution elements of the 
phase modifying element, and 



30 ^ = 0-^, and 
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K=l-J 0 {l 22x7}) , wherein 

Jo Is the zero-order Bessel function, and 

rj relates the radius Ri of the zero-order filtering region to the radius R 2 of the main- 
lobe of the Airy function of the input aperture, tj = R t 1 1^ = (0.61)" 1 ArA/, , 

5 selecting, for each resolution element, one of two phasor values which represent a 
particular grey level, and 

supplying the selected phasor values to the respective resolution elements 
(x, y) of the first phase modifying element, and 

supplying selected phasor values e^'^to respective resolution elements (x\ y') of a 
10 second phase modifying element having a plurality of Individual resolution elements 
(x\ /), each resolution element (x\ y') modulating the phase of electromagnetic 
radiation incident upon it with the respective phasor value e W Y) for generation of the 
output field o(x\yy^. 

The mathematical expressions will be further explained below. 

1 5 The axis of propagation of a plane electromagnetic field is perpendicular to the 
electric and magnetic fields. 

It should be noted that, in each resolution element of the first phase modifying 
element, one of two phasor values which represent a particular grey level of the 
amplitude component of the electromagnetic field ofr'y) may be selected. 

20 In an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields Incident upon it outside the phase shifting 
regions, i. e. A is equal to one or approximately equal to one. 

In an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields incident upon it inside the phase shifting 
25 region, I.e. B is equal to one or approximately equal to one. 

It is also preferred that the phase shift value 9 substantially fulfils the equation 

1 



K\a\ = 



2|3in|| 



for a lossless filter with A = 1 and B = L In a preferred embodiment of the present 
invention, the phase shift 9 is equal to n or approximately equal to n. Accordingly the 
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previous equation leads to Ka = / 2 and the phase values, H x > >0 . of the first phase 
modifying element may be calculated in accordance with 

where A is the Illuminated area of the first phase modifying element. 

5 The electromagnetic field or radiation may be of any frequency range of the 

electromagnetic spectrum, i e the gamma frequency range, the ultraviolet range, the 
visible range, the infrared range, the far infrared range, the X-ray range, the 
microwave range, the HF (high frequency) range, etc. The present invention is also 
applicable to particle radiation, such as electron radiation, neutron radiation, etc. 

1 0 Preferably, the electromagnetic fields are monochromatic or quasi-monochromatic so 
that the energy of the electromagnetic fields is concentrated in a narrow frequency 
bandwidth. Since the phase contrast generated amplitude pattern is reconstructed by 
interference of two electromagnetic fields generated by different phase shifting of 
different parts of the incoming field, it Is required that the frequency range of the 

1 5 emitted electromagnetic field is sufficiently narrow to ensure that the two 

electromagnetic fields are coherent so that their superposition generates the desired 
amplitude pattern, if the frequency range is too broad, the two fields will be 
incoherent and the phase information will be lost as superposition of non-coherent 
fields results in a summation of the intensities of the two fields. It is required that the 

20 difference between individual delays of electromagnetic fields to be superpositioned 
is less than the wavelength of the fields. This is a relaxed requirement that allows the 
electromagnetic fields to be relatively broad-banded. For example in the visible 
range, a Xe-lamp or a Hg-lamp can be used as a light source in a system according 
to the present invention with the advantage compared to a laser light source that 

25 speckle noise is reduced. The requirements of the spatial coherence of the 
electromagnetic fields depend upon the space bandwidth product of the 
corresponding system and how close the required system performance is to the 
theoretically obtainable performance of the system. 

Preferably, the electromagnetic radiation is generated by a coherent source of 
30 electromagnetic radiation, such as a laser, a semi-conductor laser, a strained multi- 
quantum well laser, a vertical cavity surface emitting laser (VCSEL), a maser, a 
phase-locked laser diode array, a light emitting diode, a pulsed laser, such as a sub- 
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picosecond laser, etc, or an array of such sources. However, as already mentioned, 
a high-pressure arc lamp, such as an Hg tamp, a Xe lamp, etc, may also be used and 
even an incandescent lamp may be used as a source of electromagnetic radiation. 

Each phase modifying element changes the phase of an electromagnetic field 
5 incident upon it. Optionally, it may also change the amplitude of an electromagnetic 
field incident upon it. Each phase modifying element may transmit or reflect the 
incident electromagnetic field. Each phase modifying element may be divided into a 
number of resolution elements, each of which modulates the incident electromagnetic 
field by changing its phase by a specific predetermined value. The predetermined 
10 values are assigned to each resolution element In different ways depending upon the 
technology applied in the component For example in spatial light modulators, each 
resolution element may be addressed either optically or electrically. The electrical 
addressing technique resembles the addressing technique of solid-state memories in 
that each resolution element can be addressed through electronic circuitry to receive 
15 a control signal corresponding to the phase change to be generated by the 

addressed resolution element. The optical addressing technique addresses each 
resolution element by pointing a light beam on it, the intensity of the light beam 
corresponding to the phase change to be generated by the resolution element 
illuminated by the light beam. 

20 Spatial phase modulation may be realized utilizing a fixed phase mask, a liquid 
crystal device based on liquid crystal display technology, a MEMS (micro electro- 
mechanical system), a MOEMS (micro opto-electro-mechanical system), such as a 
dynamic mirror device, a digital micro-mirror array, a deformable mirror device, etc, a 
membrane spatial light modulator, a laser diode array (integrated light source and 

25 phase modulator), smart pixel arrays, etc. 

Seiko-Epson produces a transmitting liquid crystal SLM (LC-SLM) having a high 
resolution matrix of transparent liquid crystal elements wherein the relative 
permittivity of each element can be electrically modulated in order to vary the 
refractive index and thereby the optical path length of the element, 

30 Meadowlark produces a parallel-aligned liquid crystal (PAL-SLM) with a high fill 
factor, but this device has a very low resolution in that it contains only 137 phase- 
modulating elements. 

Hamamatsu Photonics produces a dynamically controllable PAL-SLM with VGA or 
XGA resolution. 
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Texas Instruments produces a Digital Mirror Device (DMD) having an array of mirrors 
each of which can be tilted between two positions. 

The spatla) phase fitter is typically a fixed phase mask, such as an optically flat glass 
plate coated with a dielectric layer In the region wherein the modulated 
5 electromagnetic field is phase shifted Q in relation to the remaining part of the 
electromagnetic field. However, the spatial phase modulators mentioned in the 
previous section may also be used for spatial phase filters. In addition, non-linear 
materials providing self-phase modulation, such as Kerr-typa materials, can also be 
used for introducing the phase shift 6. 

10 An imaging system maps the phase modulating resolution elements of the first phase 
modifying element onto the second phase modifying element This imaging system 
may comprise a 4Mens configuration (two Fourier transforming lenses utilizing 
transmission of light or one Fourier transforming lens utilizing reflection of light) or a 
single imaging lens. However, any optical imaging system providing a filtering plane 

15 for the spatial phase filter may be applied In a phase contrast imaging system. 

In the method and system according to the present invention, the electromagnetic 
field o(x\ y 1 ) is generated by superposition of electromagnetic fields in the image 
plane of the imaging system, The first phase modifying element changes the phase 
values of an electromagnetic field incident upon it and the imaging system directs the 

20 electromagnetic field with changed phases reflected from or transmitted through the 
phase modifying element towards the spatial phase filter. The phase filter phase 
shifts a part of the electromagnetic field and the Imaging system is adapted to 
superimpose in the image plane the phase shifted part of the electromagnetic field 
with the part of the electromagnetic field that is not phase shifted by the spatial phase 

25 filter. 

According to a preferred embodiment of the invention, the first phase modifying 
element is positioned at the front focal plane of a lens white the spatial phase filter is 
positioned in the back focal plane of the lens, whereby a first electromagnetic field at 
the phase modifying element is Fourier transformed by the lens into a second 

30 electromagnetic field at the phase filter. Thus, specific spatial frequencies of the first 
electromagnetic field will be transmitted through the spatial phase filter at specific 
positions of the phase filter. For instance, the energy of the electromagnetic field at 
zero frequency (DC) is transmitted through the phase filter at the Intersecting point of 
the Fourier plane and the optica! axis of the lens also denoted the zero-order 

35 diffraction region. 
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it is an advantage of the invention that utilisation of arrays of sources is facilitated in 
that the positioning and/or shapes of the phase shifting regions of the phase filter 
may be matched to the geometry of the source. For example, if a linear array of 
VCSEls forms the source, the phase shifting regions of the spatial phase filter form a 
5 corresponding linear array of phase shifting regions, each of the regions being 
positioned at the zero-order diffraction region of a respective VCSEL in the VCSEL 
array. Further, the shape of each phase shifting region may match the shape of the 
zero-order diffraction region of the respective VCSEL 

Likewise, a phase filter may match a source with a specific geometrical shape with a 
1 0 continuous phase shifting region covering an area of the phase filter that corresponds 
to the zero-order diffraction region of the source. 

Thus, the energy of the electromagnetic fields of the system may be distributed over 
a large area compared to the area of a zero-order diffraction region of a single plane 
electromagnetic field of a known phase contrast imaging system. 

1 5 Thus, the phase shifting regions of the spatial phase filter may form a rectangular 
array, a circular array, a linear array, two linear crossing arrays, a continuous region, 
a ring, eta 

At least two substantially plane electromagnetic fields with different axes of 
propagation may be generated in a time multiplexed manner, e.g. by a scanning 
20 mirror or prism, deflecting or reflecting a beam of electromagnetic field in different 
directions of propagation. 

The capability of handling high energy levels of electromagnetic fields of the present 
invention may be utilized for provision of a 3D laser cutter comprising a system 
according to the present invention. 

25 Further, the capability of handling high energy levels in combination with the 
capability of generating a desired three-dimensional field comprising desired light 
beams may be utilized for provision of an optical tweezer or an array of optical 
tweezers according to the present invention. 

In an embodiment of the present invention, wherein the apertures of the system is of 
30 insignificant importance to the operation of the system and calculation of the phasor 
values, K is equal to one or approximately equal to one. 
For a more detailed understanding of the invention, theZernike approximation is 
reviewed below, followed by a generalization where the above-mentioned 
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mathematical expressions are derived for an on-axis centred phase contrast filtering 
Implementation. 

The Zernike phase contrast method allows for the visualization of phase 
perturbations by the use of a Fourier plane phase shifting filter. The Dutch physicist 

5 Fritz Zernike received the Nobel Prize In 1953 for inventing this method, which led to 
a break-through in medicine and biology by making essentially transparent cell or 
bacteria samples clearly visible under a microscope. Its successful operation, 
however, requires that the spatial phase distribution, #{x,y) , at the input is limited to 
a "small-scale" phase approximation where the largest phase is typically taken to be 

10 significantly less than n /3-According to this assumption, a Taylor expansion to first 
order is sufficient for the mathematical treatment so that the input wavefront can be 
written as 

*xp{?$(x*y))*i+%(x,y) O) 

The light corresponding to the two terms in this "small-scale" phase approximation 
15 can be separated spatially by use of a single lens where the phase distribution is 
located in the front focal plane and the corresponding spatial Fourier transformation 
is generated in the back focal plane of the lens. In this first order approximation the 
constant term represents the amplitude of on-axis light focused by the lens in the 
back focal plane and the second spatially varying term represents the off-axis light 
20 Zernike realized that a small phase shifting quarter wave plate acting on the focused 
light makes it possible to obtain an approximately linear visualization of small phase 
structures by generating interference between the two-phase quadrature terms in Eq. 
(1): 

/(x' >: >/'H + 2*(x',;k') (2) 

25 It should be noted that a three-quarter waveplate works equally well to produce 
contrast, but the plus sign in Eq. (2) is negated leading to so-called negative phase 
contrast. A substantial improvement in the visibility of the Zernike phase contrast 
visualization in Eq. (2) requires strong damping of the focused light in addition to the 
phase shift required to generate the contrast 

30 In the general case r where we are not limited to a small-scale input phase 
perturbation we cannot assume that a series expansion to first order as in the 
Zernike approximation is a sufficient representation of a given phase perturbation. 
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Higher order terms in the expansion need to be taken into account, so the expansion 
takes the form: 

(3) 

However, here the spatially varying terms can not be considered as separate from 
5 the supposedly focused light represented by the first term in this Taylor series 
expansion, as is implied by the Zemike approach, and all of these spatially varying 
terms contribute to the intensity of the on-axls focused light. For a significant 
modulation in the input phase, this contribution of the spatially varying terms can 
result in a significant modulation of the focal spot amplitude in the back focal plane of 
10 the lens. These terms can in fact result in either constructive or destructive 

interference with the orvaxis light, although the net result will be an attenuation of the 
focused light amplitude, which only has a maximum value for a perfect unperturbed 
plane wave at the input. 

For phase objects not fulfilling the Zemike approximation we must, therefore, find an 
1 5 alternative mathematical approach to that of the Taylor expansion given in Eq . (3)- 
We have chosen a Fourier analysis as a more suitable technique for completely 
separating the on-axis and higher spatial frequency components. This gives the 
following form for exp(/$*,y)), where (x,y) eH: 

exp(/<j>(*,.y))- l\dxdy\ JJexp^*,;/))^ (4> 
\n J q 

20 In this Fourier decomposition the first term is a complex valued constant linked to the 
on-axis focused light from a phase object defined within the spatial region, Q , and 
the second term describes light scattered by spatially varying structures in the phase 
object. Comparing Eq. (3) and Eq. (4) it is apparent that the first term of Eq. (3) is a 
poor approximation to the first term of Eq. (4) when operating beyond the Zemike 

25 small-scale phase regime. 

An important Issue to consider when analysing the effect of spatial filtering of the light 
diffracted by phase perturbations is the definition of what spatially constitutes focused 
and scattered light. In the previous description of Zemike phase contrast it was 
assumed that the focused light is spatially confined to a somewhat unphysical delta 
30 function. As we know, any aperture truncation inherent in any practical optical system 
will lead to a corresponding spatial broadening of the focused light. It is therefore 
essential that we define the terms "focused llghf and "scattered light" explicitly for 
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such a system. In this context it is necessary to look more carefully at the sequence 
of apertures confining the light wave propagation through a typical optical set-up. 

A commonly applied architecture that provides an efficient platform for spatial filtering 
is illustrated in Pig. 1 and is based on the so-called 4-f configuration. An output 

5 interferogram of an unknown phase object or phase disturbance is obtained by 
applying a truncated on-axis filtering operation In the spatial frequency domain 
between two Fourier transforming lenses (L1 and L2). The first lens performs a 
spatial Fourier transform so that directly propagated light is focused Into the on-axis 
filtering region whereas spatially varying object information generates light scattered 

1 0 to locations outside this central region. We can describe a general Fourier filter in 
which different phase shifts and amplitude damping factors are applied to the 
"focused" and "scattered 0 light In Fig 1 , we show a circularly symmetric Fourier filter 
described by the amplitude transmission factors A and 8 for the "scattered 11 and 
"focused" light respectively and by the relative phase shift 6, These filter parameters 

15 can be chosen to replicate any one of a large number of commonly used filter types 
(i.e. phase contrast, dark central ground, point diffraction and field-absorption 
filtering). By applying a given Fourier fitter and a second Fourier lens, we obtain an 
interference pattern in the observation plane. The focused on-axis light acts as the 
synthetic reference wave (SRW) in the common path interferometer (CPI) system, 

20 this interferes with the scattered light to generate the output interference pattern. In 
the following section we discuss the importance of the SRW and show how it 
influences, among other things, the choice of the Fourier filter parameters. 

Having described the generic optical system that makes up the CPI, we turn to a 
detailed analytical treatment of the important elements in this system. Assuming a 
25 circular Input aperture with radius, Ar , truncating the phase disturbance modulated 
onto a collimated, unit amplitude, monochromatic field of wavelength, X , we can 
describe the incoming light amplitude a(x f y) by, 

a (xy) =circ&/Ar ) exp (fy (x, y)) (5) 

at the entrance plane of the optical system shown in Fig. 1 using the definition that 
30 the clrc-function is unity within the region, r * ffi+y* £ Ar , and zero elsewhere. 
Similarly, we assume a circular on-axis centred spatial filter of the form. 



H{f x , /,) ~a[\ + (BA-* exp (/9)-l)circ(/,/A/;)] 



(6) 
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where B e=[0; l] is the chosen filter transmittance of the focused light, 0 e[0;2;r] is 
the applied phase shift to the focused light and A e[0;l] is a filter parameter 
describing field transmittance for off-axis scattered light as indicated in Fig, 1 . The 
spatial frequency coordinates are related to spatial coordinates in the filter plane 

5 such that: (/„/,) = (#T'(W/) and X = Jfi + ft ■ 

Performing an optical Fourier transform of the input field from Eq. (5) followed by a 
multiplication with the filter parameters in Eq. (6) and a second optical Fourier 
transform (corresponding to an inverse Fourier transform with inverted coordinates) 

we obtain an expression for the intensity /(x'.y^K'^Of describln 9 the 
10 interferogram at the observation plane of the 4-f set-up: 

I{x\f)=\ \A 2 \exp{iH*'>y))^ < 7 > 

where g(r') is the synthetic reference wave (SRW) and the terms a and ?(x' f y) 
are given by: 

15 (8) 

It should be noted that to achieve a tractable analytic expression In Eq. (7) it has 
been assumed that the spatial frequency content of the phase object is sufficiently 
described by the term, a , within the on-axis centred filtering region characterized by 
20 the spatial frequency range A/ f - 

The generally complex valued and object dependent term, a , corresponding to the 
amplitude of the focused light plays a significant role in the expression for the 
interference pattern described by Eq, (7). Referring to the discussion in the 
introduction, we are now able to confirm that the frequent assumption, that the 
25 amplitude of the focused light is approximately equal to the first term of the Taylor 
expansion in Eq. (1), can generally result in misleading Interpretations of the 
interferograms generated at the CPI output. 

Of similar importance in the analysis of Eq. (7) is the term g(r') describing the spatial 
profile of the SRW, diffracted from the aperture formed by the on-axis centred filtering 
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region. It is the Interference between this SRW term, carrying the information about 
the filtering parameters, and the imaged phase object that generates the output 
interferogram. Thus, it is important to obtain an accurate description for the SRW and 
thereby an accurate derivation for Eq. (7). The zero-order Hankel transform followed 
5 by a series expansion in the spatial dimension, f, will be used to describe the SRW. 
This is a relatively simple approach, which to the knowledge of the author has not 
previously been applied to this problem. 

For a circular input aperture with radius, Ar , we can describe the radius of the 
corresponding central phase shifting region of the Fourier filter (characterized by the 
1 0 parameters 3 and 0) in terms of a radial spatial frequency range A/, . We can thus 
obtain the following expression for the SRW by use of the zero-order Hankel 
transform: 

g(r <) = 2*Ar f J x (2*Ar/ r }/ 0 (2v'/r VK (9) 

In order to simplify the analysis, we introduce a term ij, which explicitly relates the 
15 radius of the central filtering region, fl, , to the radius of the main-lobe of the Airy 
function, R x , resulting from the Fourier transform of the circular input aperture alone. 
We can thus express tj in terms of Ar and tf r such that: 

T\ = RjRi={p6\ttotf, (10) 

where the factor of 0.61 arises from the radial distance to the first zero crossing of 
20 the Airy function corresponding to half of the Airy mainlobe factor, of 1 .22. If we make 
this substitution in Eq.(9) and then perform a series expansion in r' , we obtain the 
following expression for the SRW: 

g(rO-l-y o (l.22in,)-[(0.61nn)V 2 (1.22nn)](rVAr) J + 

{[(0.6l7Oi) 3 /4][2/ 3 (]^2im)-0,61nn/ 4 0 - 22l "l)]} 0* *V ° ^ 

25 

In this expansion, the SRW is expressed in radial coordinates normalised to the 
radius of the imaged input aperture. This can easily be scaled to allow for a 
magnification within the imaging system, though for the remainder of our analysis a 
direct imaging operation is assumed. From Eq.{1 1) it is apparent that the SRW will 
30 change as a function of the radius of the central filtering region. Additionally, it is 
clear that the SRW profile is not necessarily flat over the system output aperture, 
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Thls is an important, yet often neglected factor in determining the performance of a 
CPI. 

Depending on the accuracy needed for the description of the interferograms one can 
choose to include a number of spatial higher order terms from the expansion in Eq. 
5 {1 1). The influence of the higher order terms has the largest impact along the 
boundaries of the imaged aperture. For 77 -values smaller than 0.627 and when 
operating within the central region of the image plane, spatial higher order terms are 
of much less significance and we can safely approximate the synthetic reference 
wave with the first and space invariant term in Eq. (1 1): 

10 g(r<<=central region) *\~J 0 {L22nx\) (12) 

so that we can simplify Eq- (7) to give: 

/(x^')s^|exp(/$(^\yO)+^R(^ ,e ^( /9 )- 1 )f < 13 > 

where K = 1 - J q {122tt7j) . The influence of the finite on-axis filtering radius on the 
focused light is thus effectively included as an extra "filtering parameter* so that the 
15 four-parameter filter set {A.BAKfy)) together with the complex object dependent 

term, a , effectively defines the type of filtering scheme we are applying. 

Having determined a suitable operating range for the CPI in terms of the production 
of a good SRW, we must now examine the role that the remaining filter parameters 
play in the optimisation of a CPI. 

20 From Eq, (13) we see that the filter parameters {A.B.O) can be combined to form a 
single complex valued term, C r the combined fitter term, such that: 

C = [C|exp(/M/ C ) = BA' X exp (i9)-1 (14) 

therefore, Eq. (13) can be simplified to give: 

J(^yW a |exp^ (15) 

25 where 



BAT 1 - ^l+2|C|cos{^ c )+|C| l 
sme*(M- l )~ X \C\$\n(yr c ) 



(16) 
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Since it is a complex variable, the combined filter term C, which effectively describes 
the complex filter space, can be considered to consist of a vector of phase y/ c and 
length |C| as shown in Eq.(14). Thus in order to obtain an overview of the operating 
space covered by all the possible combinations of three independent filter 
5 parameters (A, B,0) we can now instead choose to consider a given filter In terms of 
the two combined parameters y/ c and |c| . However, referring to Eq.(15), it can be 
seen that the filter parameter, A also appears independently of the combined filter 
term, C. Fortunately, this issue can be resolved by considering that the term BA" } 
from Eq. (14) must be constrained in the following way; 

10 

BA~ X < 1 => A-\ B=\C+\\ 

BA"=\ => ,4 = 1, £ = 1 (17) 
M ml >\ => j*»l f A~\C+\[ l 

These constraints arise from the adoption of a maximum irradlance criterion 
minimising unnecessary absorption of light in the Fourier filter, which reduces both 
1 5 irradiance and the signal to noise ratio in the CPI output 

In the previous sections we derived expressions relating the spatial average value of 
a given phase disturbance to obtain peak irradlance and optimal visibility in 
combination with high accuracy in systems with unknown wavefront phase 
disturbances. We saw that if a CPI is applied to wavefront sensing or the 

20 visualisation of unknown phase objects the Generalised Phase Contrast (GPC) 
method specifies the filter phase and aperture size parameters for achieving optimal 
performance in extracting and displaying the phase information carried by the 
incoming wavefront On the other hand, in cases where we have control over the 
incoming wavefront or phase modulation the GPC method provides extra means of 

25 optimisation by encoding the phase distribution itself In addition to modifying the filter 
parameters. The two main scenarios: A) synthesizing the spatial phase for intensity 
display or B) measuring the spatial phase with high accuracy, strongly influences 
which of the parameters In the analysis that should be kept fixed and which could be 
changed or adapted. The first approach is particularly useful when the filter 

30 parameters have a restricted dynamic range or are fixed. The rigorous derivation of 
the equations for choosing these parameters will be derived In this section. 

When synthesizing an input phase distribution for optima! visibility of an output 
intensity pattern the situation Is more relaxed than the situation involving accurate 
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interferometric measurements of unknown phase disturbances. The parameter 77 can 
therefore in most cases be chosen to completely encompass the zero-order light with 
the result that the term, K , tends to unity as the Bessel function tends to zero In Eq. 
(12). For this particular case, the SRW becomes a flat top profile where we can 
5 achieve nearly 100% light efficiency. For smaller and irregular phase patterns fine- 
tuning of 7 in the region 0.4-0.6 provides for an efficient operating regime while 
maintaining minimal losses. 

In order to optimize a synthesized light distribution for maximum contrast, we wish to 
generate an intensity distribution with a lowest intensity equal to zero, i.e. in at least 
10 one point (xj,j>j): 

/(<.jfcft)-0 (20) 

where fa is the relative phase shift generating zero-intensity in (xj 9 yf) of the 
observation plane. 

Applying this dark background condition in Eq. (13) we can obtain the following 
15 expression for a no-loss phase-only filter with filter transmission 
parameters,^ 5 = 1: 

^P(l-exp(ie))=exp(/| 0 ) (21) 

A key point arising from Eq. (21) is that we now have a simple way of expressing a 
new design criterion relating the spatial average value of any input phase pattern to 
20 the zero-order phase shift of a matched Fourier phase filter. Since K is by definition 
positive and by taking the modulus of Eq.. (21) we obtain: 

iCpl^^sinte/a)!" 1 (22) 

Eq. (22) is a key result for the fully transmlssive wavefront engineered GPC mapping 
that makes it possible to deduce the range of valid phase parameters fulfilling our 

25 design criteria from Eq (20). The largest possible value that the term, JC|a|, takes on 

is unity, this leads to the following solution interval for Eq. (22) within a full phase- 
cycle: 

9 = [rc/3;5ji/3] (23) 



From Eq. (22) we also observe that /den can take on a value limited to the interval: 
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*H-l* i] (24) 

Eq. (22) and the solution intervals described by Eqs. (23)-(24) specify the design 
parameters for achieving optimal performance in extracting and displaying the phase 
information carried by the incoming wavefront. Moreover, Eq. (22) hints towards extra 
5 means of optimisation by encoding the phase modulation depth itself in addition to 
the no-loss phase-only filter. This last approach is particularly useful when the filter 
phase has a restricted dynamic range or is fixed. 

Now, assuming that we have a fixed and fully transmissive phase-only filter, the best 
choice for the filter parameter is a value that allows for the largest dynamic range of 
1 0 phasor values at the input. Accordingly, the smallest possible ma! value, K<x - 1 / 2 , 
is desirable implying that 0 = k , leading to the output intensity distribution: 

/(^y)=2[l-co 5 (4»(^,/))] (25) 

Inserting Ka - 1 / 2 in Eq. (8) we obtain the following two requirements for the input 
encoded phase function $(x t y) : 



15 



A (26) 
KA~ l l\sut(<l>(x,y))dxdy=0 



where A is the illuminated area of the phase modifying element 

We observe that it is only the first requirement in Eq. (26) that is directly related to the 
output intensity in Eq. (25) via the cosine term. Since there are always two choices 
for a given phasor value that result in the same cosine value (excluding 0 and n) t we 
20 notice that the second requirement can subsequently be fulfilled independently of the 
first requirement simply by complex conjugating an appropriate number of phasor 
values. This fact Is a key feature of the GPC-method since it makes it possible to 
solely concentrate on the first requirement in the process of synthesizing a desired 
an virtually no-loss grey level intensity pattern. 

25 The first requirement in Eq. (26) can be fulfilled by several means, including, dynamic 
phase range adjustment, fill factor encoding, phase-histogram adjustment, spatial 
scaling of phasor pattern, raster encoding etc. In a histogram adjustment technique 

one will typically start out with a desired relative intensity distribution /(x'.y)^,^ 

where the maximum achievable intensity level is unknown but relative Intensity levels 
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are known and the lowest intensity level is fixed by the background criterion of Eq. 

(20). The procedure is now to adjust the histogram for /(V,/)^^ while 

maintaining identical relative intensity level ratios until the first requirement in Eq. (26) 
is fulfilled Subsequently, the second requirement in Eq. (26) is fulfilled by complex 

5 conjugating an appropriate part of the phasors. The simplest procedure is to complex 
conjugate every second identical phasor value independently of the spatial location. 
However, this "phasor flipping** procedure can also be turned into an advantageous 
tool (an extra degree of freedom) for manipulating the spatial frequency content in 
order to optimize the separation of low and high spatial frequency terms at the 

1 0 Fourier filter plane by taking the spatial phasor location into account E.g. 

neighbouring phasor values can be chosen to have a maximum difference between 
them, thereby Introducing high spatial frequency modulation easing the filtering in the 
spatial Fourier domain. 

In most cases, however, equalized output intensity levels are sufficient in the 
15 succeeding analysis, we therefore focus on the encoding of the input phase levels to 
achieve binary output intensity levels. A derivation based on ternary phase levels 
allows for the widest range of binary intensity pattern encoding and automatically 
provides for the simplified but important binary phase level encoding as a special 
case. 

20 For the ternary phase encoding, we consider the illuminated portion of the input 
aperture area, A % as divided into sub-areas A 0l A t and A 2 with respective phase 
values $ 0 , $ x and 0 2 . We are aiming for the derivation of general expressions 
relating the addressing parameters for the phase modulation to the range of possible 
phase parameters of the Fourier filter obeying the design criterion we have already 

25 set out* We can express the total truncated area and its average phase modulation, 
as the sum of the phase-weighted sub-areas: 

A 0 exp(^)+ A, exp(/^)+ A 2 exp(rft) - A|a| (27) 

This can be further simplified by expressing the sub-areas as fractions of the total 
area, A, such that F t = A r / A and F l « A 2 / A : 

30 (l-F } -FjJexp^oJ-f/JJexp^J+^exp^^pl (28) 

As previously mentioned we are interested in binary intensity patterns with levels 
corresponding to the Input phase values. In this case the dark background region Is 
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10 



15 



defined by (a 0 ,^) and the bright output level of intensity, /, is determined by 

(A„4) and (A 2> &) at the Input plane. For me binary output intensity condition it 
follows that 



/($,) = /(*,) 



(29) 



5 This equality corresponds to a symmetric condition that can be easily verified 
applying the phasor chart analysis technique demonstrated in section 4. Due to this 
symmetry we can simplify the analysis by applying the following substitution: 

A* = $,-$ 0 =4o-*2 (30) 
so that Eq.(28) can be rewritten as: 

F^exp^Hj+F,^ (31) 

It is now a straightforward task to solve Eq- (31) for the real part and the imaginary 
part respectively, to obtain the following set of equations: 



F x + F 2 = (2it - l)(2is: (l - cos ( A4>))) M 
F x - F 2 = sin (8)(2tf sin (A*)(l - cos (D)))" 



This can also be expressed in terms of the fractional areas, such that: 



F x * (4K)- ] [(2/T - - cos^))" 1 +sin (9) (sin (A*) (l -cos(9)))~ ! ] 



(32) 



*i=(4*)" 



(2AT - l)(l -cos(A4)f ' - sin (9)(sm (A*)(l - cos(e)))"' ] 



(33) 



20 Since we have focused on solutions where identical intensity levels are obtained in 
both the F r region and the F r region we can define the resulting illumination 
compression factor, C, in the following way: 



C^(F^F 2 y l -(l-(2i:r , )" I (l-cos(A(j))) 



(34) 



The minimum compression factor corresponds to uniform illumination at the output 
25 such that FJ + F 2 = 1 , whereas the maximum compression factor is found to be 
C-»oofor tf=^. 
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An Interesting special case can be deduced from Eq. (31) by setting F 2 m 0, where 
we find that: 

F-Ft=[K(l-w{»))-l)(^^ (35) 
Implying that for the binary phase modulation case we must have: 

5 A<j> = 9 (36) 

In order for the fill factor, F, to be real-valued. 

This result turns out to be the special case that corresponds to the set of solutions 
where a binary phase pattern serves as the input. 

The second phase modifying element imposes a spatial phase shift, \y(x\y'), on the 
1 0 phase contrast generated electromagnetic field o(x\ y 1 ) . This superposition generates 
an arbitrary controllable complex field, with amplitude given by the 

square-rooted Intensity, (*',/), and phasor given by the remaining 

part under division with the amplitude term, o(x , 9 y f )e ili '' y) /\o(x l i y')\ m The resulting 

arbitrary controllable complex field, o (x^y)*** 1 ^ is capable of re-distributing the 
1 5 light into an arbitrary three-dimensional focusing within a selected volume of 

operation. This focusing within a volume Is a result of complex wave propagation and 
can be deduced by use of Maxwell's equations. A simplified scalar description of this 
complex wave propagation into a 3D field distribution can be obtained by use of a 
simple plane wave Fourier decomposition: 

operating on the Fourier transform, 3 , of the controllable complex electromagnetic 
field leaving the second phase modifying element Any subsequent focusing optics 
can be included in the Fourier decomposition of Eq. (37). 

At least one of the first and second phase modifying elements may further be 
25 adapted for phase modulation by first phasor values for a first polarization and 
second phasor values for a second orthogonal polarization of the Input 
electromagnetic field. Individual phase modulation of orthogonal polarizations of an 
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etectromagnetic field may for example be performed by a birefringent spatial light 
modulator, such as a spatial light modulator based on liquid crystal technology- 

This allows for an extended functionality where the generated 3D field distribution 
can be divided into two orthogonal and non-interfering polarisation components that 
5 can e.g. be applied in a counter-propagating geometry 

Preferably, the second phase modifying element is adapted for phase modulation by 
first phasor values e 1 * 1 **-* for a first polarization and second phasor values e m ^ tV) 
for a second orthogonal polarization of the input electromagnetic field. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 Fig. 1 illustrates schematically a known 4f phase contrast imaging system, 

Fig. 2 illustrates schematically an embodiment of the present invention, and 

Fig. 3 illustrates schematically an embodiment of the present invention generating 
the desired three-dimensional field by focusing of o(x\ y')^*-* 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

15 Fig- 1 shows a known 4f CP! phase contrast imaging system 1 . A laser (not shown) 
emits a light beam, which is expanded by a beam expander (not shown) into a plane 
light wave of uniform intensity and directs it towards a first phase modifying element 
4. The light beam Is transmitted through the first phase modifying element 4 and a 
Fourier transforming lens LI. The first phase modifying element 4 is positioned in the 

20 front focal plane of the lens L1 and a spatial phase filter 6 is positioned in the back 
focal plane of the lens L1 that is also the front focal plane of a lens L2. The Fourier 
transforming lenses L1 , L2 need not have identical focal lengths. Different focal 
lengths lead to a magnification ratio different from one. The spatial phase filter 6 
phase shifts by ft and optionally attenuates (by a factor B), the zero order diffraction 

25 part 8 of the light phase modulated by the first phase modifying element Optionally, 
the remaining diffraction part of the light modulated by the phase modifying element 
may be attenuated by a factor A. The electromagnetic field o(x\ /) is generated in 
the back focal plane 9 of the lens L2. 

Figs, 2 and 3 illustrate schematically an embodiment to of the present invention, 
30 comprising a 4f CPI phase contrast imaging configuration as illustrated In Fig. 1. 
Corresponding parts in the figures are designated with identical reference numerals. 
It is obvious for the person skilled in the art that the 4f configuration may be 
substituted by the 2f or 1f configurations disclosed in WO 96/34207. The operation of 
the 4f CPI phase contrast imaging system is explained with reference to Fig. 1 and is 
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not repeated.. Again the electromagnetic field o(x\ y') Is generated In the back focal 
plane 9 of the lens L2 In front of the second phase modifying element 14. The first 
phase modifying element 4 has a plurality of individual resolution elements (x, y), 
each resolution element (x, y) modulating the phase of electromagnetic radiation 
5 incident upon it with a predetermined phasor value e^K As already disclosed, the 
phasor values ef 0 *** of the phase modifying element and the phase shift value 9 are 
calculated substantially in accordance with 



10 A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction region, 

B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region, 



a a ja|exp(/^) is the average of the phasors e Wx,y) of the resolution elements of the 
1 5 phase modifying element, and 



JC=sL-J 0 (l .22^) , wherein 

J 0 is the zero-order Bessel function, and 

tj relates the radius R» of the zero-order filtering region to the radius R 2 of the main- 
20 lobe of the Airy function of the Input aperture, 1] = $/^= (0.61)" 1 brbf f . 

For each resolution element (x, y) t one phasor value is selected from two phasor 
values representing a particular grey level. 

The computer 12 supplies the selected phasor values e^ to the respective 
resolution elements (x, y) of the first phase modifying element 4, and supplies the 
25 determined 6 value to the spatial phase filter 6. 

The second phase modifying element 14 also has a plurality of individual resolution 
elements (x\ y'), each resolution element (x*, y') modulating the phase of 
electromagnetic radiation incident upon it with a predetermined phasor value e^*'^. 
The computer 12 supplies determined phasor values e^^to respective resolution 
30 elements {x\ y') of the second phase modifying element 14 for modulation of the 




wherein 



^ = and 
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phase of the electromagnetic field o{x\ y l ) incident upon it with the respective phasor 
value 8^ y| for generation of the desired output field o(x\ y^e 1 *** 0 ^ emitted from the 
second phase modifying element 14. 

Thus, any desired amplitude as a function of (x\ y') and any desired phase as a 
5 function of {x\ y") of the output field o(x\ y'Je 1 ^- ^ may be synthesized by the 
method and apparatus according to the present invention. 

The output field o(x', y , )e w ' t *'* yj propagates and generates the desired three- 
dimensional field u(x", y", z"), Le. the light is re-distributed into an arbitrary three- 
dimensional field distribution within a selected volume. A simplified scalar description 
10 of this complex wave propagation into a three-dimensional field distribution can be 
obtained by use of a simple plane wave Fourier decomposition: 



where 3 is the Fourier transformation of the controllable complex field o(x\ y*)e w <*fl 
leaving the second phase modifying element 14. 

1 5 A further optical system, such as a lens 16, a microscope objective lens, a curved 
mirror, an aspheric lens, etc, may focus the electromagnetic field o(x\y , )e lvt * 10 within 
the selected volume. Any subsequent focusing optics can be included in the Fourier 
decomposition of the above equation. 

Further, the computer 12 may comprise light control means for controlling the power 
20 of the light-emitting source generating the field incident on the first phase modifying 
element 4. The computer may also comprise input means, such as a keyboard, a 
mouse, a 3D mouse, 3D virtual reality equipment, a diskette drive, a USB interface, 
an optical disc drive, a network interface, a modem, etc, for receiving a three- 
dimensional field to be synthesized by the system 10, From the received three- 
25 dimensional field, the computer may be adapted to calculate phasor values and 
eF* rt to be transmitted to the respective resolution elements {x, y) and (x\ y') of the 
first and second phase modifying element 4, 14, and the phase shift 0 of the spatial 
phase filter 6 for transmission to the spatial phase filter 6 In accordance with the 
above-mentioned equations. 




dfjfy 



30 
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PROPOSED APPLICATIONS: 

• 3D optical micro* and nanomanipufation in reaMime. 

• 3D optical multi-beam tweezing for manipulation of micro-objects, such as micro- 
components, biological cells, etc, using electromagnetic gradient forces 

5 proportional to the optical intensity pointing in the direction of the intensity 
gradient, 

• Optical fractionation, sorting, sifting etc. 

• Efficient and dynamic spot-array generators to provide bias or holding beams for 
3D arrays of photonic elements, such as bistable elements, photonic switches 

10 and smart pixels. 

• Generation of structured light for machine vision applications. E.g. periodic and 
skew periodic mesh grid illumination in 3D that can be updated in parallel 

» Photolithographic applications (laser 3D direct writing in parallel without the need 
for sequential scanning). E.g. high power laser direct writing of waveguides in 
15 Ge-doped silica. 

• Volume light intensity modulation in general by use of pure phase modulation 
(radiation focusators). 

• 3D laser beam shaping in real time. 

• 3D image projection without the need for a laser-scanning device. 
20 • Dynamic Infrared Scene Projection (D1RSP). 

• Exposure device for grating and mask production. 

• LlDAR applications. 

• Laser printing In parallel. 

• Laser show applications. 
25 • Atmosphere research. 
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CLAIMS 

1 A phase contrast system for synthesizing an output electromagnetic field 
u(x", y", z")» comprising 

a first phase modifying element for phase modulation of an input electromagnetic 
5 field by phasor values e l ^ t 

first Fourier or Fresnel optics for Fourier or Fresnel transforming the phase 
modulated electromagnetic field positioned in the propagation path of the phase 
modulated field, 

a spatial filter for filtering the Fourier or Fresnel transformed electromagnetic 
10 radiation by 

in a region of spatial frequencies comprising DC In the Fourier or Fresno] 
plane 

phase shifting with a predetermined phase shift value 0 the modulated 
electromagnetic radiation in relation to the remaining part of the 
15 electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant 0, and 

in a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
20 with a constant A 

second Fourier or Fresnel optics for forming an electromagnetic field o(x\ y') by 
Fourier or Fresnel transforming the phase shifted Fourier or Fresnel transformed 
electromagnetic field, and 

a second phase modifying element for phase modulating the electromagnetic field 
25 o(x\ y') into the electromagnetic field o(x\ y'le 14 ^ * propagating as the desired 
output electromagnetic field . 

2. A phase contrast system according to claim 1 , wherein at least one of the first and 
second phase modifying elements is further adapted for phase modulation by first 
phasor values for a first polarization and second phasor values for a second 
30 orthogonal polarization of the input electromagnetic field. 
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3. A phase contrast system according to claim 2, wherein the second phase 
modifying element is further adapted for phase modulation by first phasor values 
e wicx\ jo f 0r a first polarization and second phasor values e ,^ ' 2^ * 4V ' , for a second 
orthogonal polarization of the input electromagnetic field. 

5 4. A phase contrast system according to claim 2 or 3, further comprising an element 
for directing the phase modified orthogonal fields into separate paths of propagation, 
e.g, to be applied in a non-interfering counter-propagating geometry. 

5. A phase contrast system according to any of the preceding claims, wherein 

A=1. 

10 6. A phase contrast system according to any of the preceding claims, wherein 

B = 1. 

7. A phase contrast system according to any of the preceding claims, wherein 

e=ir. 

8. A phase contrast system according to any of the preceding claims, wherein the 
15 phasor values e^ of the first phase modifying element and the phase shift value 9 

substantially fulfil that 

a(^/)s^[exp(^ 

wherein 

A Is an optional amplitude modulation of the spatial phase filter outside the zero- 
20 order diffraction region, 

S is an optional amplitude modulation of the spatial phase filter In the zero-order 
diffraction region, 

a = (ajexp(^) is the average of the phasors e^ of the resolution elements of the 
phase modifying element, and 
25 f$ = ^-^,and 

K = 1 - J 0 (1 .22;n?) f wherein 

J 0 is the zero-order Bessel function and 
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r\ relates the radius of the zero-order filtering region to the radius R 2 of the main- 
lobe of the Airy function of the input aperture, r\ = //^ » (0.6l)~ l ArA£ . 

9. A phase contrast system according to any of the preceding claims, wherein the 
phase shift value 9 substantially fulfills the equation 

1 0. A phase contrast system according to any of the preceding claims, wherein at 
least one of the first and second phase modifying element comprises a complex 
spatial electromagnetic field modulator that is positioned in the path of the Input 
electromagnetic field and comprises modulator resolution elements (Xm, y m ), each 

10 modulator resolution element (x ro , y m ) modulating the phase and the amplitude of the 
electromagnetic field incident upon it with a predetermined complex value 

1 1 . A phase contrast system according to any of the preceding claims, further 
comprising a light source for emission of the input electromagnetic field, the light 

15 source comprising a laser array, such as a VCSEL array. 

12. An optical micro-manipulation or multi-beam optical tweezer system according to 
any of the preceding claims. 

13. A laser machining tool according to any of the preceding claims. 

14. A method of synthesizing an output electromagnetic field u(x n , y'\ z"), comprising 
20 the steps of 

phase modulating an input electromagnetic Held by phasor values e^* ,v> , 

Fourier or Fresnel transforming the phase modulated electromagnetic field, 

filtering the Fourier or Fresnel transformed electromagnetic radiation by 

In a region of spatial frequencies comprising DC In the Fourier or Fresnel 
25 plane 

phase shifting with a predetermined phase shift value 0 the modulated 
electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
30 with a constant B, and 
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in a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constants, 

forming an electromagnetic field o(x\ y 1 ) by Fourier or Fresnel transforming the phase 
5 shifted Fourier or Fresnel transformed electromagnetic field, and 

phase modulating the electromagnetic field o(x\ y') into the output electromagnetic 
field o(x\ y'Je^-fl > propagating as the desired output electromagnetic field 
u(x n ,y n ,z u ). 

15. A method according to claim 14, further comprising the steps of 

1 0 dividing the electromagnetic field ofx*,/) into pixels in accordance with the disposition 
of resolution elements (x, y) of a first phase modifying element having 



1 5 calculating the phasor values e^ of the phase modifying element and the phase 
shift value 9 substantially In accordance with 



wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
20 order diffraction region, 

B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region. 

a = ja|exp(/^) is the average of the phasors e 1 ^ of the resolution elements of the 
phase modifying element, and 
25 <£=<j>-<jh,and 



a plurality of individual resolution elements (x, y), each resolution element 
(x, y) modulating the phase of electromagnetic radiation incident upon it with 
a predetermined phasor value a****. 




K = 1 - J 0 (1 ,22n7]) , wherein 

Jo is the zero-order Bessel function, and 
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ij relates the radius Ri of the zero-order filtering region to the radius R2 of the main- 
lobe of the Airy function of the input aperture,^ = ^ IR^ « (O^l)" 1 ArA£ , 

selecting, for each resolution element, one of two phasor values which represent a 
particular grey level, and 

5 supplying the selected phasor values e* M to the respective resolution elements 
(x, y) of the first phase modifying element, and 

supplying selected phasor values e^'^to respective resolution elements (x\ y') of a 
second phase modifying element having a plurality of individual resolution elements 
(x\ /), each resolution element (x\ y*) modulating the phase of electromagnetic 
10 radiation incident upon it with the respective phasor value e w * ,y) for generation of the 
output field o(x\y')e w ** 
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Fig. 1 
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Fig. 2 
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